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Abstract. In the present work, we investigate the formation of Stokes profiles 
and spectro-polarimetric diagnostics in an active region plage near the limb. We 
use 3-D radiation-MHD simulations with unipolar fields of an average strength 
of 400G, which is largely concentrated in flux tubes in which the field reaches 
typical kilo-Gauss values. We generate synthetic Stokes spectra by radiative 
transfer calculations, then we degrade the simulated Stokes signal to account for 
observational conditions. The synthetic data treated in this manner are com- 
pared with and found to roughly reproduce spectro-polarimetric high-resolution 
observations at u=0.39 obtained by the SOUP instrument with the Swedish 1-m 
Solar Telescope at the beginning of 2006. 



1. Introduction 



The study of small scale magnetic flux concentrations near the solar limb pro- 
vides additional insight in to the existing models of magnetic flux concentra- 
tions ( Frutiger et al. 20031 ). These flux concentrations have been associated 



with solar faculae. Understanding the physical processes be hind solar faculae 



through observations and simulations i s an important topic (ILites et al. 200 



Keller et al. 200i [Carlsson et al. 20041: ISteiner 20051 : Iffirzberger k Wiehr 200_. 
Okunev Kneer 20051 : iDe Pontieu et al. 2006l h One of the motivations for 
studying the facular phen omenon is the influen ce of their brightness on the 
Sun's irradiance variation ( Fligge Solanki 200 ll ). 



In order to gain further insight into magnetic flux concentrations in plage 
regions at intermediate [i values, we perform a study of the spectropolarimet- 
ric signal in MHD simulations and observations originating from plage regions 
at n = 0.39. The simulations used here ar e fully compressible 3-D radiation- 
MHD simulations ( Vogler Sz Schiissler 20031 : IVogler et al. 20051). The observa - 



tions has been recorded with the SOUP instrument ( Title Rosenberg 1 98 ll ) 
at the Swedish 1-m Solar Telescope (SST) in 2006. 
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Figure 1. Magnetic field strength map at an optical depth T5000 = !• 



2. Methods and simulations 



The 3D radiation-MHD simulations used include, a solution of the fully com- 
pressible MHD equations, including parti al ionization in the equat i on of state 

and a non-local non-grey radiative transfer ( Vogler Schiissler 20031 ; IVogler et al. 20051 ). 
The MHD simulation box is 6x6 Mm in the horizontal plane and 1.4 Mm in the 
vertical direction, with a resolution of 288x288x100 grid points. The mean mag- 
netic field strength is about 400 G. (See Figures [U and EJ). In these Figures we 
can see the fully developed granulation pattern which has interacted with the 
magnetic field. The flux gets concentrated in intergranular lanes as a result of 
flux expulsion ( Schiissler 19901 ). The flux concentrations appear darker than the 
average granular intensity. This is due to the partially suppressed plasma mo- 
tion through field lines. Nevertheless, a proper analysis of energy exchange via 
convection and radiation is necessary before dra wing conclusions on the thermal, 
and radiative properties of flux concentrations ( Schiissler et al. 2 003). This is 
beyond the scope of the present paper. 

We construct an inclined view of the simulation box by interpolation of the 
different quantities along rays with the desired inclination angle. This allows 
comparisons with observations away from disk center. Fig. [3] shows an inclined 
view of Fig. [2] at an angle of 67 degrees (/i = 0.39). One can notice the higher 
values of intensity at some locations (e.g. at coordi nates (2200km, 30 0km) and 
(3000km, 600km)). These are identified as faculae dKeller et al. 20041 s ). 

In order to compare these 3D radiative-MHD simulations with Spectropo- 
larimetric observations, we have to calculate the Stokes signal emerging from 
the inclined simulation box (e.g. the continuum intensity map near 6302.5 A 
is shown in Figure [3]). This is done using the STOPRO code in the SPINOR 
package (Solanki 1987; iFrutiger et al. 2000 ). In order to take into account ef- 
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Figure 2. Normalized continuum intensity map near 6302.5 A emerging 
from the simulation box at disc center. 
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Figure 3. Normalized continuum intensity map near 6302.5 A emerging 
from the simulation box at \x = 0.39. 



fects such as the finite aperture of the telescope, and seeing, we perform then a 
spectral and spatial smearing of these synthetic polarimetric data. 



3. Comparison with spectropolarimetric observations at Mu=0.39 

We compare the emerging Stokes signal from the simulation run with a set of 
high spatial resolution spectropolarimetric images obtained with the SOUP in- 
strument at the Swedish 1-m Solar Telescope. The observational data consists 
of 9 images: the four Stokes polarization images at both + and -50 mA from the 
line center of the Fe I 6302.5 A line, and a wide band image. The images are 
restored using the Multi-O bject Multi-Frame Blin d Deconvolution (MOMFBD) 
image restoration method ( van Noort et al. 20051 ). The SOUP images are based 



on 480 exposures at each line position which means that the effective exposure 
time amounts to 7.2 s. The combined use of the SST adaptive optics system 
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and MOMFBD post-processing resulted in a spatial resolution in the Stokes ob- 
servations that approaches the diffraction limit of the telescope (better than 0.2 
arcsec). In order to find the actual resolution of the polarimetric observations, 
we carry out the following steps: 1/ We choose a quiet-sun simulation run and 
convolve the synthetic Stokes profiles in the spectral dimension with a Lyot type 
filter of FWHM=70 mA (similar to SOUP). 2/ In order to account for the lower 
spatial resolution in the observations, we apply a low-pass filter to the synthetic 
images which has the shape of a top-hat function and effectively removes power 
at the highest spatial frequencies - beyond the spatial resolution of the observa- 
tions. In addition, we convolve the synthetic images with a Lorentzian profile 
which accounts for the far wings of the PSF that are not corrected for in the 
MOMFBD restoration. 3/ We pick a quiet region from the observed Stokes-/ 
images either at + 50 mA or -50 mA and compare it with the synthetic ones. 
The matching of the two data sets is done by comparing their power spectra 
(Figure H] upper panel) and their standard deviations. This is done through 
an iterative process where both the FWHM of the PSF and Lorenz profile are 
allowed to change until the two data-sets fit with each other. 
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Figure 4. Two lower panels: Stokes-I maps at +50 mA from line center Fe 
I 6302.5 A. Lower left: simulations, lower right : observations. The upper 
panel shows the integrated power spectra for the two lower images (full line : 
observations, dashed line: simulations). 



The so obtained degrading parameters are used when comparing more active 
regions, like the ones in Figure The two lower panels of Figure 0] indicate the 
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Figure 5. Upper left panel: observed Stokes- V at -50 mA from line center 
Fe I 6302.5 A. Lower left panel: Stokes-V from simulations degraded to the 
resolution of observations. A measure of the Stokes-V signal at the white 
slits-positions is shown on the right panels. 



observed (right one) and degraded-simulated (left one) Stokes-I maps at +50 
mA from line center (Fe I 6302.5 A). The simulated map corresponds to a quiet 
region where the mean field strength is about 20G. One can see the similarity of 
the two images. The original observed image was inclined by an angle of about 
20 degrees, which we have corrected here. 

In the two left images in Figure we see an observations and a degraded 
simulation of Stokes-V at -50 mA. A direct quantitative comparison between 
Stokes-V signals in these two images can be done by plotting the signal along 
slits (e.g. the white slits in the left panels). The resulting signal is shown in the 
two right panels. We observe that the two curves are reasonably similar. This 
indicates that the observed features here are well reproduced by an MHD simu- 
lation run with an average field strength of 400G. At this stage we can address 
the question: does the MHD simulation-run reproduce the actual photospheric 
situation at the observed region on figure [5]?. The results above are in favor of 
a positive answer. But the fact that we have only two spectral positions at +50 
mA and -50 mA makes it necessary to proceed to further comparisons with a 
larger spectral sampling. Then, in case simulations and observations match, we 
can proceed further by studying these regions in detail directly from the MHD 
simulations. 
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4. Conclusions 

We presented here a comparative study about small-scale magnetic flux con- 
centrations near the solar limb. This is part of a more general work aiming at 
investigating their physical and spectropolrimetric properties through Stokes- 
diagnostics and direct MHD simulations. We focussed here on identifying sim- 
ilarities between observed and simulated Stokes signal at + and -50 mA from 
line center of Fe I 6302.5 A. The comparison has been made possible, by in- 
troducing suitable instrumental and seeing degradation to the simulated Stokes 
signal. We thus identified similarities between the observed and simulated maps 
of Stokes-V at two wavelength positions + and -50 mA from line center of Fe 
I 6302.5 A. Although we see similarities, we mention the necessity of further 
investigation before concluding that the simulations are indeed reproducing the 
observed magnetic features. Such investigations should include more spectral 
positions, allowing a more consistent study of the variation of physical quanti- 
ties (e.g. temperature, magnetic field vector, line-of-sight velocity) through the 
inclined magnetic flux concentrations. 
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